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ABSTRACT 

We present an analysis of the multiwavelength behaviour of the blazar OJ 248 at 
z = 0.939 in the period 2006-2013. We use low-energy data (optical, near-infrared, 
and radio) obtained by 21 observatories participating in the GLAST-AGILE Support 
Program (GASP) of the Whole Earth Blazar Telescope (WEBT), as well as data from 
the Swift (optical-UV and X-rays) and Fermi (y-rays) satellites, to study flux and 
spectral variability and correlations among emissions in different bands. We take into 
account the effect of absorption by the damped Lyman a intervening system at z = 
0.525. Two major outbursts were observed in 2006-2007 and in 2012-2013 at optical 
and near-IR wavelengths, while in the high-frequency radio light curves prominent 
radio outbursts are visible peaking at the end of 2010 and beginning of 2013, revealing 
a complex radio-optical correlation. Cross-correlation analysis suggests a delay of the 
optical variations after the y-ray ones of about a month, which is a peculiar behaviour 
in blazars. We also analyse optical polarimetric and spectroscopic data. The average 
polarization percentage P is less than 3%, but it reaches ~ 19% during the early stage 
of the 2012-2013 outburst. A vague correlation of P with brightness is observed. There 
is no preferred electric vector polarisation angle and during the outburst the linear 
polarization vector shows wide rotations in both directions, suggesting a complex 
behaviour/structure of the jet and possible turbulence. The analysis of 140 optical 
spectra acquired at the Steward Observatory reveals a strong Mg II broad emission 
line with an essentially stable flux of 6.2 x 10” 15 ergcm~ 2 s -1 and a full width at 
half-maximum of 2053 km s _1 . 

Key words: galaxies: active - galaxies: quasars: general - galaxies: quasars: individ¬ 
ual: OJ 248 - galaxies: jets. 
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Figure 1. Light curve of OJ 248 in R band. It includes 1211 points from 15 observatories, distinguished by different colours and symbols. 
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1 INTRODUCTION 

The emission of active galactic nuclei (AGNs) is explained 
by the existence of a supermassive black hole (SMBH) at 
the centre of the host galaxy, which converts gravitational 
energy of material located in the surroundings into electro¬ 
magnetic energy. This material forms a disk and loses angu¬ 
lar momentum due to the viscosity in the disk, finally falling 
onto the black hole. In general, the AGN spectra may show 
broad and narrow emission lines produced in regions close 
to the nucleus. Sometimes they can also show lines from the 
host galaxy. In radio-loud AGNs two plasma jets are ejected 
in direction perpendicular to the disc. 

Among the different types of radio-loud AGNs, the ob¬ 
jects called “blazars” (BL Lacs and flat spectrum radio 
quasars, FSRQs) are powerful emitters from radio wave¬ 
lengths up to 7 -ray energies. They present strong flu x var i- 
ability and high and variable polarization fe.g. ISmi SdH). 
The most accepted scenario to explain these features sug¬ 
gests that we are observing the emission from a jet of ma¬ 
terial accelerated to relativistic velocities in the vicinity of 
the SMBH, and oriented very close to our line of sight. Thus 
the jet radiation is Doppler boosted and dominates over the 
other emission components from the nucleus (disc, broad 
line region - BLR, narrow line region) or host galaxy. The 
origin of the low-frequency radiation (radio to UV or X- 
ray band) from the jet is attributed to synchrotron emission 
and the high-energy radiation (X- to 7 -rays) to an inverse- 
Compton process by the same relativistic electrons pro¬ 
ducing the synchrotron photons. After the launch of satel¬ 
lites for high-energy observations such as the Astrorivelatore 
Gamma ad Immaqini Leqqero ( AGILE: iT avani e t al.ll200fll j 
and Fermi dAbdo et al.ll2009l : I Atwood et al.ll2009l b the num¬ 
ber of sources detected at 7 -rays has increased significantly, 


allowing a more detailed investigation of the high-energy 
processes occurring in blazars. 

Among blazars, in general BL Lacs are objects of lower 
luminosity with featureless spectra or very weak emission 
lines. In contrast, FSRQs have higher luminosities and 
stronger emission lines. In this paper, we present multifre¬ 
quency observations of the FSRQ OJ 248 (0827+243) in 
2006-2013 performed in the framework of a campaign led by 
the Whole Earth Blazar Telescope^] (WEBT). The WEBT 
radio-optical observations are complemented by high-energy 
data from the Swift and Fermi satellites. 

In t he Roma BZ CAT multi-frequency catalog of 
blazarfl dMassaro et al.l [20091 ) OJ 248 appears with a red- 
shift a = 0.939 flagged as uncertain. Mg II a nd Fe II absorp¬ 
tion l ines at z = 0.525 were detected by lUlrich &; Owenl 
Jl977il in the source optical spectrum. This intervening 
Damp ed Lyman a. (DLA) sys tem was subsequently stud¬ 
ied bv lRao fc Turnshekl (2000), who estimated a hydrogen 
column density IVh = (2.0 ± 0.2) x 10 20 cm -2 . The DLA 
syste m is likely a disc galaxy with on going star forma¬ 
tion (jSteidel el al.lljQO^ : lllao et al.ll2003l b Although it does 
not affect the b lazar photometry because of its faintness 
dRao et al.ll2003h . its absorption of the source radiation must 
be taken into account. 

OJ 248 was detected by the Energetic Gamma Ray 
Experiment Telescope (EGRET) instrument on board the 
Compton. Gamma Ray Observatory (CGRO) with a vari¬ 
able flu x. In the third EGRE T catalog of high-energy 7 -ray 
sources dHartman et al.|[l999h it appears with a flux F(E > 
100 MeV) = (24.9+3.9) x 10 -8 ph cm -2 s _1 , single measure- 

1 http://www.oato.inaf.it/blazars/webt 

2 Edition 4.1.1, August 2012; http://www.asdc.asi.it/bzcat 
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Figure 2. Light curves of OJ 248 in optical (R band) and near-IR 
(J, H , K bands) in 2006-2013 built with GASP-WEBT data. In 
the near-IR light curves the blue diamonds are from the Campo 
Imperatore Observatory and the red points from the Teide Ob¬ 
servatory. 


merits ranging from 15.6 to 111.0 x 10 -8 phcm _ 2 s _1 . It has 
F(E > 100 MeV) = (5.3 d= 0.5) x 10 -8 phcm -2 s _1 and a 
spectral in dex 2.67 d= 0.07 in the Second Fermi-LAT catalog 
(2FGL; iNolan et al.ll20l3) . In the Nine ties the source was 
very active also in the optical band (e.g. IVillata et al.lfl997l : 
IRaiteri et alJll998h . 


2 OPTICAL AND NEAR-IR PHOTOMETRY 

Optical photometric observations were provided in the R 
band by several observatories participating in this WEBT 
project, including the GASP-WEBT collaboration and the 
Steward Observatory progra m in support of the Fermi 7 -ray 
telescope rtSmith et al. 120091 1. They are: Abastumani (Geor¬ 
gia, FSU), Calar Altc^ (Spain), Crimean (Ukraine), Low¬ 
ell (Perkins telescope, USA), Lulin (Taiwan), Mt. Maidanak 
(Uzbekistan), Roque de los Muchachos (Liverpool telescope, 
Spain), Rozhen (Bulgaria), San Pedro Martir (Mexico), St. 
Petersburg (Russia), Steward (USA), Teide (IAC80 tele¬ 
scope, Spain), Tijarafe (Spain), Torino (Italy), Vidojevica 
(Serbia). 

The period of interest goes from 2006 March up to 2013 
July. We collected a total of 1356 data points, 1211 of which 

3 Calar Alto data was acquired as part of the MAPCAT project: 
http: //www. iaa.es/~iagudo/research/MAPCAT 


survived the light curve cleaning process, through which we 
discarded data with large errors as well as clear outliers. 

In the light curve of Fig. [T] we can see two major flaring 
periods, in 2006-2007 and 2012-2013. The first flare is of 
similar b rightness (R = 15.6) as that observed in November 
1995 by IRaiteri et al.1 (|l998f ). which was a historical max¬ 
imum. The second outburst appears more prominent and 
shows a stronger variability. Two minor events are visible in 
2009-2010 and in early 2011. 

The 2007 peak was very sharp and was characterized 
by a brightening of ~ 1.3 mag in 11 days and about 1.0 mag 
fading in 9 days. A noticeable variability also characterises 
the 2012-2013 outburst, with variations of about 1 mag in 


f 

Raiteri et al. 

il998|), of 

Villata et al. 

19971), and 

Fan et al. (2004)'). But the 


I. 16 mag in 63 days observed by 
of 1.05 mag in 58 days reported by 
source also exhibits short-term variability. I 11 particular, we 
found a couple of changes of ~ 0.3 mag in about 7.5 hours 
in late 2012. I ntrad ay variability was previously reported by 
IRaiteri et al.l d 1998lj . who observed a brightness decrease of 
0.73 mag in 20 hours. 

The GASP-WEBT near-IR data are collected in the 

J, Ft, I\ bands at the Campo Imperatore (Italy) and Teide 
(TCS, Spain) observatories. Deta ils on the da t a acq uisi- 
tion and reduction are given in IRaiteri et al.l (12014 1. A 
comparison between the R-band and near-IR source be¬ 
haviour is shown in Fig. [5] Although the near-IR light curves 
are less sampled than the optical one, we can recognize 
the same main features, in particular the two outbursts 
of 2006-2007 and 2012-2013. The variation amplitude in- 
crea ses with wavelength , as it is usually observed in FSRQs 
(e.g. IRaiteri et al]|2012l l. suggesting the presence of a “sta¬ 
ble” blue emission component, likely thermal radiation from 
the accretion disc. 


3 OPTICAL POLARIMETRY 

Blazars are known to show variable polarization in both po¬ 
larized flux percentage (R) and electric vector polarization 
angle (EVPA) (e.g. ISmithlfl996l J . In particular, wide rota¬ 
tions of the linea r polarisation vector have been detected in a 
number of cases (Marschcr ct al. 2008, 2010; Lar ionov et al.l 
l2013l : ISorcia et akll2014l l. which have been interpreted as the 
effect of motion along spiral trajectories. 

Optical polarization data for this paper were provided 
by the Calar Alto, Crimean, Lowell, San Pedro Martir, St. 
Petersburg, and Steward observatories. In Fig. [3] we show 
the time evolution of the polarization percentage P com¬ 
pared with the R-band light curve. For most of the time, the 
source showed low P (average value of ~ 3%), but during 
the brightening phase of the 2012-2013 outburst P reached 
~ 19%, suggesting a correlation between P and brig htness 
typical of FSRQs (Ag. lSmithlll996l : IRaiteri et al .1120 1 .4) . Fig¬ 
ure U shows P versus the de-absorbed flux density in the R 
band (see Sect. 8 ). For any value of Fr there is a large disper¬ 
sion of P, but the highest values of P (> 13%) are reached 
when the source is bright (Fr > 1.3mjy). The linear Pear¬ 
son’s correlation coefficient is 0.60, indicating a marginal 
correlation. 

In order to investigate whether possible rotations of the 
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Figure 3. Time evolution of the optical magnitude in R band 
(top panel) and of the percentage of polarized flux (bottom 
panel). The data are from different observatories: Calar Alto 
(green squares), Crimean (red diamonds), Lowell (pink asterisks), 
San Pedro Martir (cyan triangles), Steward (blue circles), and St. 
Petersburg (orange crosses). 


linear polarisation vector occurred, we first examine the be¬ 
haviour of the Q and U Stokes’ parameters. Figure [o] shows 
that OJ 248 spent most of the time with Q and U being close 
to zero. With this being the case, even small variations in 
Q and U necessarily lead to large EVPA rotations that are 
difficult to accurately follow unless the observations are very 
dense. The problem is mitigated during the 2012-2013 out¬ 
burst, when Q and U exhibit large variations. This appears 
clearer in the Q versus U plot in Fig. [6] where subsequent 
data belonging to short time periods with good sampling 
have been connected to show the time evolution during the 
2012-2013 outburst. With this in mind, in Fig. [7] we finally 
plot the EVPA as a function of time during the outburst. 
The ± nn ambiguity was fixed by assuming that the most 
likely value is that minimizing the angle variation, i.e. we 
added/subtracted 180° when needed to minimize the dif¬ 
ference between subsequent points separated by less than 5 
days. It seems that there is not a preferable EVPA value and 
that the linear polarisation vector underwent wide rotations 
in both directions and with different ranges of angles. All 
this suggests a complex behaviour/structure of the magnetic 
field in the jet as is expected e.g. from a turb ulent plasm a 
flowing at a relativistic speed down the jet IIMarscheiil20 14). 


4 OPTICAL SPECTROSCOPY 

Several previous studies found that there is no correla¬ 
tion between the jet activity and the b ehaviou r of the 
broad emission line s in blazars (see e.g. ICorbett et al.l[2000l : 
iRaiteri et al.ll2007l ). Indeed, the BLR gas is likely ionized 
by the accretion disc radiation, whose variability is weaker 
and occurs on longer time sca les than the varia bility of the 
beamed jet emission (see e.g. iKaspi et al.l l2000). However, 
iLeon-Tavares et all d2013h detected a flare-like variability of 
the Mg II emission line in the blazar 3C 454.3 during an 



Figure 4. Optical polarization percentage P plotted against the 
de-absorbed flux density in the R band. The data are from dif¬ 
ferent observatories: Calar Alto (green), Crimean (red), Lowell 
(pink), San Pedro Martir (cyan), Steward (blue), and St. Peters¬ 
burg (orange). 
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Figure 5. The evolution of the Q and U Stokes’parameters as a 
function of time. Different colours during the 2012-2013 outburst 
highlight the data of the selected periods indicated in Fig. [6] 


outburst and claimed that the broad emission line fluctua¬ 
tions are linked to the non-thermal continuum emission from 
the jet. With this in mind, we analyse the spectroscopic be¬ 
haviour of OJ 248 during our monitoring period. 

Optical spectra were taken at the Steward Observatory 
of the University of Arizona for the “Ground-based Obser¬ 
vational Support of the Fermi Gamma-ray Space Telescope” 
prograirQ. Data for this program are tak en at the 2.3 m Bok 
telescope and 1.54 m Kuiper telescope llSmith et al.l 120091 ); 
140 spectra of OJ 248 were acquired during the first five cy¬ 
cles of the Fermi mission, from October 2008 to June 2013. 

All spectra show a prominent Mg II AA2796, 2803 broad 
emission line, which was measured by fitting a Gaussian 
model with a single component, after subtracting a linear 
fit to the continuum. We used ad-hoc routines developed in 


4 http://james.as.arizona.edu/~psmith/Fermi/ 
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Figure 6. Q versus U for all the data shown in Fig.[5| Coloured 
lines connect subsequent data belonging to the periods listed in 
the legend (JD-2450000). The direction is indicated by the arrows. 
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Figure 7. Time evolution of the optical magnitude in R band 
(top panel), of the percentage of polarized flux (middle panel) 
and of the EVPA “corrected” for the d=n7r ambiguity (bottom 
panel, see text for explanation) during the 2012-2013 outburst. 
Data are from different observatories: Calar Alto (green squares), 
Crimean (red diamonds), Lowell (pink asterisks), San Pedro Mar- 
tir (cyan triangles), Steward (blue circles), and St. Petersburg 
(orange crosses). 


IDL and based on the MPFH0 libraries llMarkwardtjl2009l 'l . 
The continuum region was selected from adjacent regions to 
the emission line, which are free of features. The uncertainty 
of the measured flux was estimated from the average of the 
residuals obtained after the continuum fitting and was then 
used to determine the parameter confidence limits applying 
the routines in the MPFIT library. 

In Fig. [ 8 ] we show two of the spectra corresponding to 
different brightness states and in Fig. [9] the line and contin¬ 
uum fits as well as the residuals. The results of the Mg II 
analysis are shown in Fig. 1101 The line flux presents some dis¬ 
persion around a mean value of 6 . 2 x 10 -15 erg cm ~ 2 s _1 with 
standard deviation of 0.5 x 10 ~ 15 erg cm -2 s -1 . The possible 
presence of line variability can be checked e.g. by calculating 
the mean fractional variation / = %/ a 2 — 8 2 /< F >, where a 
is the standard deviation, S the mean sq uare un certa inty of 
the fluxes, and < F > the average flux (|Petersonll200lh . The 
result is / = 0.08, which means that the line flux is basically 
stable, and this is true also during the 2012-2013 outburst 
period, when the continuum flux increased by a factor ~ 6 . 
Moreover, no delayed line flux increase was detected also 
after the outburst. Hence, the enhanced jet activity respon¬ 
sible for the outburst does not affect the BLR. 

We measured the line full width at half-maximurrQ 
(FWHM), from which one can derive the velocity of the 
gas clouds in the BLR. The average value is vfwhm = 
2053kms _1 with a standard deviation of ~ 310kms -1 . 
The corresponding de-projected gas velocity, of course, de- 
pends on th e geometry a nd orientation of the BLR (see e.g. 
iDecarli, Dotti fe Trevesl201llh In the blazar model, the BLR 
should be nearly face on, so measurements of FWHM are 
likely to be underestimated, since the measurement of the 
radial velocity will likely miss most of the orbital component. 

In the bottom panel of Fig.[l0]we finally plot the equiv¬ 
alent width (EW) versus the continuum flux density. The 
EW decreases when the source brightens, which confirms 
that the jet is not the ionizing source of the BLR . We no¬ 
t ice t hat, according to the classical definition dStickel et al.l 
EA blazars with rest-frame EW less than 5 A are classi¬ 
fied as BL Lacs; in the case of OJ 248, this happens when 
the observed EW goes below 5 Ax(l + z) ~ 9.7 A, which 
occurs when the source continuum flux density around the 
Mg II lin<D exceeds ~ 0.6 x 10 -15 erg cm -2 s -1 A -1 . This 
underlines the limit of the classical distinction between BL 
Lacs and FSRQs based on the EW, which depends on the 
source brightness. 

Fig. [8]shows a strong absorption line at A = 4270 A due 
to Mg II absorption in the intervening DLA system at 2 = 
0.525 mentioned in the Introduction. We cannot perform a 
more detailed analysis of this line because the resolution of 
our spectra is about 20 A, implying a velocity resolution of 
~ 1400kms _1 , while the FWHM of the Mg II absorption 
line is ~ 270kms _1 (ISteidel et al.ll2002h . 


5 http://www.physics.wisc.edu/~craigm/idl/fitting.html 

6 Corrected for the instrumental broadening of the line. 

7 We estimated the continuum flux density in the spectral regions 
5320—5360 A and 5500—5530 A and then took the mean value. 
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Figure 8. Optical spectra of OJ 248 during different brightness 
states obtained at the Steward Observatory, showing the Mg II 
emission line at z = 0.939 (blue) and the Mg II absorption line 
at z = 0.525 (orange). 
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Figure 9. Gaussian fit (red solid line) to the Mg II emission line 
of the spectrum taken on 2008 October 30 and shown in Fig. [8] 
The cyan line indicates the linear fit to the continuum; the red 
dashed line represents the residuals. 


5 OBSERVATIONS AT RADIO AND 

MILLIMETRE WAVELENGTHS 

Radio and mm observations were performed with the Medic- 
ina (5, 8, and 22 GHz), Metsahovi (37 GHz), Noto (43 GHz), 
IRAM (86 and 230 GHz), and Submillimeter Array (SMA, 
230 and 345 GHz) telescopes. 

A detailed description of the 43 GHz measurements 
performed with the Noto Radiotelescope can be found in 
Lcto ct al. (j2009j). For the Medicina observations, see e.g. 
Bach ct al. (j2007n . 

The 37 GHz observations were made with the 13.7 m 
diameter Metsahovi radio telescope. The flux density scale 
is set by observations of DR 21. Sources NGC 7027, 3C 274 
and 3C 84 are used as secondary calibrators. The error es¬ 
timate in the flux density includes the contribution from 
the measurement root mean square and the uncertainty of 
the absolute calibration. A detail ed description of t he da ta 
reduction and analysis is given in iTerasranta et alj (|l998i ). 


2009 2010 2011 2012 2013 




0.5 1.0 1.5 2.0 

Fcont [10- 15 erg cm -2 s^ 1 A -1 ] 


Figure 10. Top: time evolution of the Mg II broad emission line 
flux; the dashed line represents the average value and the dotted- 
dashed line the standard deviation around the mean. Bottom: 
EW of Mg II versus the continuum flux density; the solid line 
displays the behaviour of the EW assuming a constant line flux 
equal to its average value. 


IRAM 30 m Telescope data were acquired as part of the 
POLAMI (Polarimetric AGN Monitoring with the IRAM 
30 m Telescope) and MAPI (Monitoring AGN with Po- 
larimetry at the IRAM 30 m Telescope) programs. Data 
re duction was perform ed following the procedures described 
in lAgudo et al] d2006l . l2Q10i ). 

Millimetre and submillimetre data were also obtained at 
the Submillimeter Array (SMA) near the summit of Mauna 
Kea (Hawaii). OJ 248 is included in an ongoing monitoring 
program at the SMA to determine the fluxes of compact 
extragalactic radio s ources that can be u sed as calibrators 
at mm wavelengths llGurwell et al.l 120071 ). OJ 248 was also 
observed as part of a dedicated program to follow sources 
on the Fermi LAT Monitored Source List (PI: A. Wehrle) in 
2009 and 2010. In the ongoing monitoring sessions, available 
potential calibrators are observed for 3 to 5 minutes, and the 
measured source signal strength calibrated against known 
standards, typically solar system objects (Titan, Uranus, 
Neptune, or Callisto). In addition, from time to time cali¬ 
brator data obtained during regular science observations are 
also used to obtain flux density measurements. Data from 
this program are updated regularly and are available at the 
SMA websit«0- 

As one can see in Fig. 1111 the mm data (230 and 86 


http: //smal. sma.hawaii.edu/callist/callist.html 


© 2015 RAS, MNRAS 000.HHT61 







































Multiwavelength behaviour of the blazar OJ 248 from radio to 7 - rays 7 


2006 


2008 


2010 


2012 


q; 


15.0 

15.5 
16.0 

16.5 
17.0 

17.5 
4 
3 

3 2 
C 1 






4 


3 
S 2 

C 1 


4 


^ Q 

o, 3 

to 2 

CO 

^ 1 


4 


S 3 

m 2 

■"t 

^ 1 


4 


^ Q 

a, 3 

^ 2 

CO 

1 


4 


^ 3 

N 2 

CM 

^ 1 

4 

^ 3 

. 2 

^ 1 

4 

3 

„ 2 

^ 1 



4000 4500 5000 5500 6000 6500 

Julian Date - 2450000 


Figure 11. Light curves of OJ 248 at different frequencies in 2006-2013. From top to bottom: JJ-band optical magnitudes (see also Fig. 
[T)l. 345 GHz data from SMA, 230 GHz data from SMA (red triangles) and IRAM (blue crosses), 86 GHz data from IRAM, 43 GHz data 
from Noto, 37 GHz data from Metsahovi, 22, 8, and 5 GHz data from Medicina. 


GHz) show two prominent outbursts of the same strength 
peaking in late 2010 and early 2013. Going to lower fre¬ 
quency (37 GHz), the second outburst becomes fainter than 
the first, and it completely disappears at 8 GHz. A com¬ 
parison with the optical light curve suggests that the 2013 
mm-radio outburst is the time delayed counterpart of the 
2012-2013 optical event, while a possible correlation be¬ 
tween the optical and radio variations in correspondence of 
the 2010-2011 mm-radio outburst is more difficult to es¬ 


tablish. Indeed, there are no visible optical flares either con¬ 
temporaneous or slightly preceding the lower-frequency out¬ 
burst, but just one event about one year before, which how¬ 
ever is more likely connected with the pre-outburst bumps 
visible at 230 and 37 GHz. In fact, there is a definite rise 
in the 230 GHz light curve that starts at essentially the 
same time as the late 2009-early 2010 optical event. Most 
likely the more or less prominent optical counterpart of the 
main mm-radio outburst remained unobserved due to the 
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2010 seasonal gap: indeed, some residual activity can be 
seen at the start of the subsequent observing season. The 
2006-2007 optical outburst might be correlated with a mi¬ 
nor radio event observed a few months later at 86 (and 
230) GHz. We used the discrete correlation functio n (DCF; 
lEdelson &: Krolil3ll988l : lHufnagel fe Bregmanll 19921) to anal¬ 
yse cross-correlations among light curves. Figure [12] shows 
the DCF between the optical de-absorbed (see Sect. 0 and 
230 GHz flux densities. The peak value is 0.8, which im¬ 
plies good correlation, and the time delay corresponding to 
the peak is 28 days, which would be the time lag of the 
mm variations after the optical ones. However, this result 
is dominated by the last outburst and the DCF run on the 
pre-outburst period gives no significant signal. 

Flares can be produced by shocks propagating down¬ 
stream the jet and/or by variations of the Doppler fac¬ 
tor 5, which depends on both the bulk Lorentz factor of 
the relativistic plasma, Tb, and the viewing angle 6,5 = 
[Tb (1 — /3 cos#)] -1 , where /3 is the velocity in units of the 
light speed. The different behaviour of OJ 248 in various 
epochs, i.e. different correlation between optical and radio 
variations, can be explained in terms of a misalignment of 
the region emitting the bulk of the optical radiation with 
res pect to the zone emitting th e bulk of radio pho tons (see 
e.g. lVillata et al .1120071 . fioQQal fbl: iRaiteri et al.ll201ll ). The ra¬ 
diation coming from the jet region with a smaller viewing 
angle will in fact be more Doppler boosted. According to 
this interpretation, in 2010-2011 the radio emitting region 
was more aligned with the line of sight than the optical zone, 
while in 2012-2013 the most external jet regions, emitting 
the low-frequency radio photons, had a larger viewing angle, 
and the strong optical outburst was the effect of a viewing 
angle smaller than ever. 

We finally notice that there is a time delay of the radio 
flux variations going toward longer wavelengths; in partic¬ 
ular, by means of the DCF we could estimate the time lag 
between the 230 and 37 GHz flux changes (see Fig. US) . The 
peak of the DCF is strongly asymmetric, so that a better es¬ 
timate of the dela^ i n this case is g iven by the centroid of the 
distribution dPeterson et al.lll99a ). which indicates a value 
of 40-50 days. Moreover, the amplitude of the flux variations 
decreases at lower frequencies (the ratio between the maxi¬ 
mum and minimum flux density is ~ 24 at 230 GHz, 7.4 at 
86 GHz, and 6.3 at 37 GHz), and the light curves become 
smoother with longer-lasting events. This is what we expect 
if the radio emission at longer wavelengths comes from more 
external (because of synchrotron self-absorption) and more 
extended regions of the jet. 

An alternative picture to explain the radio-optical vari¬ 
ability is in terms of a disturbance, e.g. shock wave, propa¬ 
gating downstream the jet where the delay from high radio 
frequencies towards low radio frequencies is naturally caused 
by opacity effects. In case the optical and mm emitting re¬ 
gions are co-spatial, the reason why the mm peak is time 
delayed may be due to a rather high lower-energy cutoff to 
the electron energy distribution in the early stages of the 
outburst. 
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Figure 12. Discrete correlation functions between the R -band 
de-absorbed and 230 GHz flux densities (blue filled circles) and 
between the 230 and 37 light curves (red empty circles). 

6 SWIFT OBSERVATIONS 

The Swift satellite (jGehrels et al.l 120041 ) carries three in¬ 
struments that work simultaneously in different frequency 
ranges: the X- Ray Telescope (XRT), observing between 0.3 
and 10 keV llBurrows et al.l [20051 ). the Ultr aviolet-Optical 
Telescope (UVOT), between 170 and 600 nm dRoming et alJ 
l2005f ). and the Burst Alert Telescope (BAT), between 14 and 
195 keV jBarthelmv et al1l2005l ). OJ 248 was observed by 
Swift 86 times between 2008 January and 2013 May. 

6.1 UVOT observations 

The UVOT telescope can ac quire data in opt ical (v, b, u) 
and UV (wl,m2,w2) bands (Poole et al.112008 ). The data 
reduction was performed with the HEASoft package ver¬ 
sion 6.13 and the Calibration Database (CALDB) 20130118 
of the NASA’s High Energy Astrophysics Science Archive 
Research Centei[f] (HEASARC). We extracted the source 
counts within a 5 arcsec radius aperture and the background 
counts from a nearby circular region with 15 arcsec radius. 
We summed multiple observations in the same filter with the 
uvotimsum task and then processed them with uvotsource. 

The resulting light curves are shown in Fig. [13] No ob¬ 
servations were available in the b band. After 2010.0 the 
difference between the maximum and minimum magnitudes 
in the different bands is 1.12 in w2, 1.05 in m2, 1.15 in wl, 
1.25 in u. It is 0.67 in v, but in this case we have only four 
points. We can see that the variability in general decreases 
when the frequency increases, extending the trend we no¬ 
ticed in Sect. [2] for the optical and near-IR light curves. 

6.2 XRT observations 

Reduction of the XRT data was performed with the HEA¬ 
Soft package version 6.13 with the calibration file 20130313. 
There are 83 observations in photon-counting (PC) mode. 
We ran the task xrtpipeline with standard screening cri¬ 
teria. For further analysis we only kept the 64 observations 

9 http://heasarc.gsfc.nasa.gov 
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Figure 13. Light curves of OJ 248 built with Swift-\JVOT data 
in optical and UV. 


with more than 50 counts. Source counts were extracted 
with the xselect task from a circular region of 30 pixel 
radius centred on the source and the background counts 
were derived from a surrounding annular region of 50 and 
70 pixel radii. No correction for pile-up was needed since the 
count rate is always lower than 0.5 counts/s. We performed 
spectral fits with the Xspec package in the 0.3-10 keV en¬ 
ergy range, using the Cash statistics because of the low 
count number. We modelled the spectra with a power law 
with photoelectric absorption, adopting a hydrogen atomic 
column density Ah = 4.6 x 10 20 cm -2 obt ained by sum¬ 
ming the Galactic value 2.6 x 10 20 cm -2 (IKalberla et all 
|2005|) to that of the in tervening DLA system at z = 0.525 
llRao fc Turnshei5l2000 ). 

The X-ray spectrum acquired on 2012 October 2 is 
shown in Fig. [14] as an example. It was best-fitted with a 
power law with photon index T = 1.49 ±0.10. 

In Fig. [15] we plotted T as a function of the flux density 
at 1 keV for the 58 observations with error less than 30% 
of the flux. The values of T range between 1.06 and 2.07, 
with an average value of 1.65 and no significant trend of 
r with flux, in contrast to the harder-when-brighter trend 
sometimes found in FSRQs I'e.g. lVercellone et al.i20icl ). The 
smaller dispersion of the data points corresponding to the 
2012-2013 outburst (standard deviation a = 0.09) with re¬ 
spect to the pre-outburst data (cr = 0.24) is likely due to 
their higher precision because of the larger count number. 
In both cases the standard deviation is less than the average 


Swift-XRT 



0.5 1 2 5 

Energy (keV) 


Figure 14. The XRT spectrum of OJ 248 on 2012 October 2. 
The best fit was obtained with a power law with fixed absorption 
given by the sum of the Galactic and intervening DLA Ah values. 
The bottom panel shows the ratio of the data to the folded model. 
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Figure 15. The X-ray photon index T as a function of the unab¬ 
sorbed flux density at 1 keV. Only data with error less than 30% 
of the flux are shown. Red points correspond to data acquired 
after JD = 2456100, i.e. during the 2012-2013 outburst. 


error (0.16 and 0.29, respectively), indicating that the data 

are consistent with a constant value. _ 

Notice that ljorstad fc MarscheJ l2004h analysed Chan¬ 
dra data using Ah = 3.62 x 10 20 cm -2 , i.e. the Galac- 
tic value accordi n g to Dickev fc Lockmanl fl990), and that 
IStroh fc Falcond ( 2012 ) found a value of Ah = (7 ± 2 ) x 
10 20 cm -2 when analysing XRT data with an absorbed 
power law with freely varying Ah- 

The resulting X-ray light curve (flux densities at 1 keV) 
is shown in Fig. 1161 where it is compared to the source be¬ 
haviour at other frequencies. 


7 FERMI OBSERVATIONS 

The Fermi satellite was launched on 2008 June 11. Its aim is 
to perform a daily mapping of the 7 -ray sources in the Uni¬ 
verse. The prim ary instrument of F ermi is the Large Area 
Telescope lLAT: [Atwood et all2009h . The energy range cov- 
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Table 1. Results of the spectral analysis of the Fermi -LAT data of OJ 248 in the 0.1 TOO GeV energy range. The reference energy Eq 
was fixed to 392.1 MeV. The fitted model was a power law with photon index T. 


Date 

-^pred 

TS 

r 

Eo.l — 100GeV 
[10” 8 phcm” 2 s” 1 ] 

2008 Aug 04 - 2009 Mar 05 

539 

125 

2.77 + 0.12 

6.2 ± 0.8 

2009 Mar 05 - 2009 Oct 15 

465 

131 

2.55 ± 0.11 

5.2 ± 0.8 

2009 Oct 15 - 2010 May 27 

406 

104 

2.60 + 0.12 

5.3 ± 0.8 

2010 May 27 - 2011 Jan 05 

466 

100 

2.83 + 0.12 

6.6 ± 0.9 

2011 Jan 06 - 2011 Aug 18 

439 

97 

2.81 ± 0.13 

6.4 + 1.0 

2011 Aug 18 - 2012 Feb 02 

283 

43 

3.13 + 0.23 

5.5 ± 1.0 

2012 Feb 02 - 2012 Jul 19 

392 

107 

2.66 + 0.13 

6.9+ 1.0 

2012 Jul 19 - 2013 Jan 03 

1975 

2106 

2.38 + 0.03 

31.5 ± 1.3 

2013 Jan 03 - 2013 Jun 20 

730 

401 

2.48 + 0.07 

12.6+ 1.1 

2013 Jun 20 - 2013 Nov 08 

313 

103 

2.56 + 0.13 

6.5 ± 1.1 


ered is approximately from 20 MeV to more than 300 GeV. 
The field of view of the LAT covers about 20% of the sky, 
and maps all the sky every three hours. 

The data in this paper were collected from 2008 August 
4 (JD = 2454683) to 2013 November 8 (JD = 2456605). 
We performed the analysis with the SCIENCETOOLS soft¬ 
ware package version v9r32p5. The data were extracted 
within a Region of Interest (ROI) of 10° radius and a max¬ 
imum zenith angle of 100° to reduce contamination from 
the Earth limb 7-rays, which are produced by cosmic rays 
interacting with the upper atmosphere. Only events be¬ 
longing to the ‘Source’ class were used. The time intervals 
when the rocking angle of the LAT was greater than 52° 
were rejected. For the spectral analysis we used the science 
tool gtlike with the response function P7REP_S0URCE_V15. 
Isotropic (iso_source_v05.txt) and Galactic diffuse emission 
(gll_iem_v05.fit) components were used to model the back¬ 
ground 111 !. 

We evaluated the significance of the 7-ray signal from 
the sources within the ROI by means of the Test Statistics 
TS = 2 (logLi — logLo), where L\ and Lq are the likelihood 
of the data given the model wi th or without the source, re¬ 
spectively iMatto^etaL 1996 ). As was done in the 2FGL 
catalog dNoIan et al. KIT for the spectral modelling of 
OJ 248 we adopted a power law, N(E) = No (E/Eo)~ r , 
where Eo = 392.1 MeV is the reference energy between 0.1 
and 100 GeV. We first ran gtlike with the DRMNFB opti¬ 
mizer, including all point sources of the catalog within 15° 
from our target, and using power-law fits to model the spec¬ 
tra of these sources. We then ran gtlike a second time 
with NEWMINUIT as optimizer, after selecting the sources with 
TS > 10 and the predicted number of counts N ple d > 3. 

The results of the analysis are reported in Table Q] for 
time bins of about six months. The average photon index 
is 2.68, and its standard deviation is 0.22. This value is es¬ 
sentially the same as that reported in the 2FGL catalog 
(2.67 ± 0.07), while an analysis over the whole period con¬ 
sidered in this paper yields T = 2.56 ± 0.03 and a flux of 
(8.6 ± 0.3) x 10” 8 phcm” 2 s” 1 . The photon index variabil¬ 
ity is dominated by errors, since the variance a 2 is smaller 
than the mean square uncertainty S 2 (see Sect. 4). In Fig. 1161 


10 http://fermi.gsfc.nasa.gov/ssc/data/access/lat/ 
BackgroundModels. html 


we can see the corresponding 7-ray light curve (red points). 
We also plotted a monthly-binned light curve (blue points) 
that includes many upper limits (cyan points) because of the 
source faintness. Finally, during the 2013 outburst we per¬ 
formed weekly bins when there was a good count number to 
detail the flux variations (green points). 


8 CORRECTION FOR GALACTIC AND DLA 

EXTINCTION 

In the previous sections we presented light curves of OJ 248 
as observed magnitudes. But the near-IR, optical, and UV 
radiation from the source suffers absorption by both the 
Galaxy dust and the dust contained in the intervening DLA 
system at 2 = 0.525 mentioned in the Introduction. This is 
a problem similar to that met when ana lyzing data fro m an¬ 
other well-known blazar, AO 0235+16 (iRaiteri et al JfjQOSh . 

We estimated the Galactic reddening in the 
UBVRIJH K optical an d near-I R bands by using the 
ICardelli. Clavton fc Mathis! d 19891) extinction law, with 
Rv = Ay/E(B — V) = 3.1, which is the standard value for 
the diffuse interstellar medium. The results are reported 
in Tabled (Column 2). An estimate of the extinction due 
to the DLA system can be calculated starting from the 
hydrogen co lumn density Nh = (2.0 ± 0.2) x 10 20 cm -2 
obtained by iRao fe Turnshekl (l2000l) and adopting a gas- 
to-dust ratio equal to the average value in the Milk y Way: 
Nh = 4.93 x 10 21 cm -2 mag” 1 x E(B — V) dDiplas fe Savage! 
119941) . This yields E(B - V) = 0.04, and assuming again 
Rv = 3.1, we obtain Ay = 0. 124 at z = 0.525. Then, by 
applying the same ICardelli. Clavton fe Mathis! (Il989l ) law 
properly blueshifted, we get the values reported in Table 
0 (Column 3). We note that the DLA system is a more 
important absorber than the Galaxy. The assumption that 
the DLA system has the same absorbing characteristics 
of the Milk y Way is justified by its being a spiral galaxy 
jRao et al.ll2003f ). The values of the total extinction that 
we must apply to correct our data for absorption in both 
the Galaxy and the DLA system are given in Column 4. 

In the case of the S’uii/f-UVOT bands, because of 
the asymmetric shape of the filter responses and of the 
bumped shape of the extinction law in the UV, we calcu¬ 
lated the absorptio n in the various bands by integrating the 
ICardelli. Clavton fe Mathisl (jl989l ) law with the filter effec- 
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Figure 16. Light curves of OJ 248 at different frequencies in 2006—2013. From top to bottom: the Fermi -LAT 0.1 100 GeV fluxes 
(10“ phcm -2 s 1 ) derived with different time bins (red symbols refer to data binned over roughly six-month time intervals, blue ones 
to monthly binned data, and green symbols to weekly-binned data in the outburst period; cyan arrows indicate upper limits); the 
1-keV Swift-XRT flux densities (//Jy); the Swift-XJVOT m2 flux densities (mjy); the GASP-WEBT R-band flux densities (mjy); the 
GASP-WEBT /4-band flux densities (mjy); the 230 GHz flux densities (Jy); the 37 GHz (orange points) and 43 GHz (green points) flux 
densities (Jy). The X-ray, UV, optical and near-IR light curves were corrected by both Galactic and DLA absorption. 


live areas (see e.g. lRaiteri et al.l201(t ). The results are shown 9 CROSS-CORRELATION BETWEEN 
in Tabled VARIABILITY AT HIGH AND LOW 

ENERGIES 
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In order to better investigate the relationship between the 
source behaviour in 7-rays and that in the optical band, we 
show in Fig.[lT]the corresponding whole light curves as well 
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Table 2. Extinctio n [mag] in the various Bessel 
llBessell. Castelli fc PleJ Il998h and Swift -UVOT bands to¬ 
ward OJ 248. Both the Galactic absorption and that by the DLA 
system at z = 0.525 are given. The value of the total extinction 
suffered by the source radiation is the sum of the two. 


Band 

Ax (Gal) 

Aa(DLA) 

Total 


Swift -UVOT bands 


w2 

0.249 

0.397 

0.646 

m2 

0.261 

0.343 

0.604 

wl 

0.211 

0.331 

0.542 

u 

0.150 

0.349 

0.499 

b 

0.125 

0.240 

0.365 

V 

0.095 

0.196 

0.291 


Bessel bands 


u 

0.142 

0.316 

0.458 

B 

0.122 

0.235 

0.357 

V 

0.093 

0.195 

0.288 

R 

0.077 

0.173 

0.250 

I 

0.055 

0.131 

0.186 

J 

0.027 

0.074 

0.101 

H 

0.017 

0.045 

0.062 

K 

0.011 

0.028 

0.039 


as an enlargement of the 2012-2013 outburst period, where 
the correlation is easier to study. The start of the 7 -ray 
outburst is not covered in the optical band because of the 
solar conjunction, but the first point after the seasonal gap 
is about 0.3 mag brighter than before, suggesting that the 
outburst has already begun. In contrast, the optical light 
curve is very well sampled in the outburst decline phase, 
where the 7 - 1 'ay curve has a worst time resolution because 
of the low flux. It seems that the 7 and optical fluxes may 
well have risen and declined together, but that the period 
of major 7 activity preceded the phase of strongest optical 
activity. 

To get a quantitative estimate of this time shift, we cal¬ 
culate the DCF (see Sect. [5} between the 7 -ray fluxes and 
the R-band flux densities corrected for the total absorption 
reported in Table dE The DCF is displayed in Fig. [TS] It 
shows a peak at a lag of r p = 28 days with DCF P = 1.3 
that indicates strong correlation with the optical variations 
following the 7 -ray ones after four weeks. The delay is 29 
days if we take the centroid instead of the peak. The fig¬ 
ure inset displays the result of 2000 Monte Carlo simula¬ 
tions according to the “flux randomization-random subse t 
selection” method dPeterson et aI.fl998l : lRaiteri et al.i2003i) . 
From these simulations it is possible to estimate the uncer¬ 
tainty on the delay. We obtained that 88 % of the realizations 
led to a centroid value between 22 and 36 days. Hence, we 
infer that the optical flux variations follow the 7 flux changes 
with a delay of 29 ± 7 days. 

However, if the 7 emission is due to inverse-Compton 
scattering of soft photons off the same electrons producing 
the optical radiation, then its variations are expected to be 
simultaneous or delayed with respect to those characteris¬ 
ing the optical radiation, as resulting from modelling non- 


11 We considered the weekly 7 fluxes during the outburst period, 
and the monthly fluxes before and after that (see Fig. im while 
we binned the optical data in seven-day bins. 
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Figure 17. Light curves in the 7 -ray and optical ( R ) bands over 
the whole period (top), and during the 2012-2013 outburst period 
indicated by the yellow stripe (bottom). 


thermal flares with shocks in a jet (e. 

z. Sikora et al. 200ll: 

Sokolov. Marscher & McHardv 20041: 
200a). This was observed in several 

Sokolov & Marscher 

blazars, in particu- 


lSchniael_et_alJJ^l^), and 3C 454.3 (iBonning et al.l 120091 ; 
IVercellone et all 20 ld : iRaiteri et all l 201 lh . In contrast. 7 
variations leading the optical o nes were ob serve d e.g. 
in the FSRQs PK S 1510-09 by lAbdo et al l ll20 10l) and 
llTAmmando et ~ahl (l201ll ). and 3C 279 by Havashida et all 
1 20121) . The latter authors found a lag of about 10 days 
and explained it by assuming that the energy density of the 
external seed photons for the inverse-Compton process de¬ 
creases faster along the jet than the energy density of the 
ma gnetic fi eld causing the synchrotron optical emission (see 
also ljaniak et al1l2012l ). Analogous interpretations may also 
hold for OJ 248. 


Alternatively, the complex optical/ 7 -ray correlation 
may be explained b y con sidering the effects of turbulence 
in the jet (iMarscheri 120141. The fluctuating magnitude and 
direction of a turbulent magnetic field affects mostly the 
synchrotron radiation, and therefore adds a component to 
the optical variability that is not present in the gamma-ray 
light curve. 

Figure [IS] also displays the DCF between the 7 -ray 
fluxes and the X-ray flux densities at 1 keV, suggesting a 
delay of the X-ray variations of about 2 months. We finally 
investigated the 7 -mm correlation, finding a strong signal 
and a ~ 70 days delay of the mm variations. 

Another interesting correlation holds between the X- 
ray and mm flux densities. Indeed, beside the 2012-2013 
outburst, the X-ray light curve also shows a peak during 
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Figure 18. Discrete correlation functions between the 7 -ray 
fluxes and the JJ-band de-absorbed flux densities (blue filled cir¬ 
cles) and between the 7 -ray fluxes and the X-ray flux densities at 
1 keV (red empty circles). The inset shows the results of Monte 
Carlo simulations of the 7 -optical correlation (see text for details). 



10 BROAD-BAND SED 

Three broad-band SEDs of OJ 248 are plotted in Fig. 1201 
They correspond to the peak of the 7 -ray emission (JD = 
2456284), to the peak of the X-ray emission (JD = 2456317), 
and to a faint post-outburst epoch (JD = 2456368). The 
SEDs are built with simultaneous near-IR, optical, UV, X- 
ray, and 7 -ray data. Because of the smoother radio variabil¬ 
ity, we gave a tolerance of a few days to the radio data. 

Emission in the optical-UV receives an important con¬ 
tribution from the accretion disc radiation, whose signa¬ 
ture is more evident in the faint, post-outburst SED. The 
concave shape of the near-IR spectrum is due to the inter¬ 
section between the disc contribution and the non-thermal 
jet emissi on. These two components have been modelled by 
IRaiteri et al.l d2014l l. who found that the OJ 248 disc is more 
luminous than a typical type 1 QSO disc. 

Notice that the faintest SED at high energies has also 
the lowest radio flux at 230 GHz, but it exceeds the fluxes 
of the other SEDs at longer radio wavelengths. 
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Figure 19. Discrete correlation functions between the X-ray and 
230 GHz flux densities (blue filled circles) and between the X-ray 
and R-band de-absorbed flux densities (red empty circles). 


the maximum of the mm light curve at the end of 2010 (see 
Fig. 1161) . The DCF in Fig. 1191 displays a strong correlation 
with no time lag, which suggests that the X-ray and mm 
radiation are produced in the same region, with the X-ray 
emission likely due to inverse-Compton on the mm photons. 
A correlation between the X-ray and mm varia bility has 
already been found, e.g. in BL Lacertae (IRaiteri et al.ll2013l) . 

In Fig. llOl we also show the DCF between the X-ray and 
optical flux densities, indicating that the optical variations 
precede those in the X-rays by about one month. 

Cross-correlations of the 7 -ray, X-ray, and R-band data 
with the 37 GHz data only led to weak and somewhat con¬ 
fused signals. This is due to the different behaviour of the 
corresponding light curves, in particular to the dominance 
of the 2010-2011 outburst with respect to the 2013 one at 
37 GHz. 


11 SUMMARY AND CONCLUSIONS 

In this paper we have presented the results of a huge mul¬ 
tiwavelength observing effort led by the GASP-WEBT Col¬ 
laboration on the blazar OJ 248. Data were collected start¬ 
ing from 2006 and up to 2013, including two optical-NIR 
outbursts in 2006-2007 and 2012-2013 and two major ra¬ 
dio outbursts in 2010-2011 and 2012-2013. The 2012-2013 
outburst was also detected at high energies by the Swift 
and Fermi satellites. The correlation between the optical 
and radio outbursts is clear in 2012-2013, while the opti¬ 
cal counterpart of the 2010-2011 radio outburst is difficult 
to identify. Something likely changed in the source in the 
period between the two outbursts, one possibility being a 
slightly better alignment of the optical emitting region with 
the line of sight with the consequent increase of the Doppler 
beaming. A strong correlation between the flux variations 
at 7 -rays and those in the optical band is found, but with 
the optical variations delayed by about one month, which is 
a peculiar behaviour already found in other blazars. Strong 
correlation with no time delay has also been found between 
the X-ray and millimetre flux changes, supporting a common 
emission region in the jet. 

We have analysed the polarimetric behaviour of the 
source. The fraction of polarised flux remained low for most 
of time but during the 2012-2013 outburst, when P reached 
~ 19%. Wide rotations of the linear polarisation vector can 
reliably be detected only during the outburst and they oc¬ 
cur in both directions, suggesting a complex behaviour of the 
magnetic field in the jet possibly due to turbulence, and/or 
a complex jet structure involving spiral paths. 

Optical spectra show Mg II lines both in absorption 
and emission. The absorption line is due to an interven¬ 
ing system at z=0.525 whose reddening effects on the NIR, 
optical, and UV source emission have been estimated and 
taken into account. The presence of the intervening sys¬ 
tem must be considered also when analysing the X-ray ra¬ 
diation. As for the Mg II emission line from the source 
BLR, we estimated a mean velocity of (2053 ± 310) kms^ 1 . 
The line flux is essentially stable around a mean value of 
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Figure 20. Spectral energy distributions of OJ 248 from the radio to the 7 -ray frequencies during three epochs characterised by different 
brightness states. 


(6.2 ± 0.5) x 10” 15 erg cm -2 s _1 also during the 2012-2013 
outburst and after, confirming that the jet emission did not 
affect the BLR, even when considering up to a few months 
of possible time delay. 
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